Nuclear magnetic relaxation and superfluid density in Fe-pnictide superconductors: 

An anisotropic ±s-wave scenario 
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We discuss the nuclear magnetic relaxation rate and the superfluid density with the use of the 
effective five-band model by Kuroki et al. [Phys. Rev. Lett. 101, 087004 (2008)] in Fe-based su- 
perconductors. We show that a fully-gapped anisotropic ±s-wave superconductivity consistently 
explains experimental observations. In our phenomenological model, the gaps are assumed to be 
anisotropic on the electron- like /3 Fermi surfaces around the M point, where the maximum of the 
anisotropic gap is about four times larger than the minimum. 



o 
o 
I 

Oh 



o 

> 

0\ 



o 

00 

o 



X 



I. INTRODUCTION 

Much attention has been focused on novel Fe- 
based superconductors since the recent discovery 
of superconductivity at the high temperature 26 
K in LaFeAsOi_2;Fa;^li Up to now, many Fe-based 
superconductors (especially iron pnictides) such as 
SmFeAsOi_a;Fa; have been found and intensively 
investigated^^'^'^i^^^'ii^i^^^^ii Experi- 
mental observations of thermodynamic quantities 
and others begin now to be reported on those 
cupcrconductor: .-^^--^^-^"-^^-^^-^^-^'^-^^-^^-^'''^^-^^-'^"-'^^-^^-'^^-^'^-^^ 
Recently, the superfluid density and the nuclear 
spin-lattice relaxation rate have been analyzed 
theoreticallyi^^i^i^^i^i^ Such observations and analyses 
are important and indispensable for elucidating super- 
conducting properties, especially for Cooper-pairing 
symmetry which we will discuss. 

One of the confused points in the experiments for 
Fe-based superconductors is that the results of the nu- 
clear magnetic relaxation rate seem inconsistent with 
the superfluid density observations. The nuclear mag- 
netic relaxation rate has the lack of the coherence 
peak below and exhibits the low temperature power- 
law behavior (l/Ti cx yS)^^^^^ rpj^^^ 
ingly the evidence of unconventional superconductiv- 
ity with line-node gaps. However, some experiments 
report that the superfluid density (i.e., penetration 
depth) does not depend on the temperature at low 
temperatures, which means that the pairing symme- 
try is fully-gapped s-wave symmetryj ^^'^'^i^^i^^i^^'^'' The 
lbs-wave pairing symmetry is theoretically proposed as 
one of the candidates for the pairing symmetry in 
Fe-pnictide superconductors3ilii^ii^iMd5^^^i9 xhe 

lbs-wave symmetry means that the symmetry of pair 
functions on each Fermi surface is s-wave and the relative 
phase between them is tt. Very recently, several theoreti- 
cal groups suggested that the ibs-wave symmetry explains 
the lack of the coherence peak and the low temperature 
power-law behavior in the nuclear magnetic relaxation 
rate, with introducing impurity scatterings Part 



of their scenarios is based on the fact that, in a ±s-wave 
phase, substantial low-energy states appear in the den- 
sity of states in the case of a unitary-limit scattering, 
while only higher-energy density of states near gap edges 
is modified when approaching to the Born limit j^^^ 



To theoretically investigate the superconductivity, 
it is necessary to consider a model for the electronic 
structure. There are many theoretical studies, espe- 
cially by band calculations, to understand the unique 
electronic and magnetic properties of those Fe-pnictide 
cupcrconductor: .'^^-^-^-^^-^'^-^^-^^-^^-^''-^^-^^-^"'^^-^^-^'^-^'^-^^-^^-^''-^^-^^ 
In addition, an effective five-band model was elaborated 
by Kuroki et al.j^ where the five bands originate 
predominantly from 3d orbitals at the Fe atomic 
site. A simpler two-band Hamiltonian was also pro- 
posed as a tractable minimal model, which reproduces 
the structure of Fermi surfaces obtained by band 
calculations i^SiiliiZ^iSiZi However, Arita et al.'^^ claimed 
that the five bands are necessary for describing correct 
band dispersions around the Fermi level. They also 
suggested that an anisotropic ±s-wave superconductivity 
is realized in a Fe-pnictide superconductori^^ 



In this paper, we investigate the nuclear spin-lattice 
relaxation rate and the superfluid density on the ba- 
sis of the realistic effective five-band model. We will 
show that an anisotropic ibs-wave pair function explains 
consistently the experimental results even in assuming a 
rather clean system. 



This paper is organized as follows. The effective five- 
band model and the pair functions are introduced in Sec. 
II. We then discuss the nuclear spin-lattice relaxation 
rate (Sec. HI), the superfluid density (Sec. IV), and the 
density of states (Sec. V). Finally, the conclusion is given 
in Sec. VI. In the appendix, we describe the derivation 
of the nuclear spin-lattice relaxation rate on the basis of 
the quasiclassical theory of superconductivity. 
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FIG. 1: (Color online) (a) Band dispersion of the effective 
five-band model and (b) Fermi surfaces with the Fermi energy 
Bp = f0.97eV. 
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FIG. 2: (Color online) Schematic figures of the pair functions 
on each Fermi surface. 



II. MODEL 

We introduce the effective five-band model proposed by 
Kuroki et al^ The tight-binding Hamiltonian is written 
as 



(2.1) 



where 



creates an electron with spin a on the fi- 



th orbital at site i, niy„ = c\^^Ci^cn and t denotes 
the hopping parameters. Here, the onsite energies are 
(ei,e2,e3,e4,e5) = (10.75, 10.96, 10.96, 11-12, 10.62)eV 
and the hopping parameters are considered up to fifth 
nearest neighbors (see a table in Ref . i42il . The band dis- 
persion of this model is shown in Fig.[IJa), and the Fermi 
surfaces are shown in Fig.[ljb). There are two hole pock- 
ets (denoted as ai, 02) centered around {kx,ky) — (0,0) 
and two electron pockets around (7r,0)(/3i) or (0,7r)(/32). 

Arita et o/4^ have performed five-band RPA calcula- 
tions for the five-band model. Their result suggests that 
the pairing function is an anisotropic is-wave symme- 



try. This pairing has isotropic s-wave pair functions on 
the Fermi surfaces ai,2, and anisotropic s-wave pair func- 
tions on Pi^2 where the maximum of the pair amplitude is 
about five times larger than the minimum. Following it, 
we assume phenomenologically the anisotropic its-wave 
pair function expressed as (see Fig. ^ 



^ai.2A,2(fc) 



Ao$ci,2,ft,2(fc)tanh(aV7'c/T- 1),(2.2) 

-fa, (2.3) 
(1 + <^0min) _^ (1 - ^/3mm) COs(2(?;.i^2) 



(2.4) 



Here, $cti(2)(^) and ^/^^^^^{k) denote the pair amplitudes 
on the Fermi surfaces ai(a2) and /3i(/32), respectively. 
Equation (|2.2p with a = 1.74 reproduces well the tem- 
perature dependence of the BCS gap. The angles 0i 
and 4>2 are measured from the (tt, 0) direction around 
{kx,ky) — (tt, 0) and (0, tt), respectively. The range of the 
gap-anisotropy parameter ^pmin is < $;3min < 1- The 
larger $/3min within this range, the weaker anisotropy. 
The sign of — $a corresponds to the relative phase of the 
pair functions between the a and /3 Fermi surfaces. If $a 
is positive (negative), the pairing is ±s-wave (s-wave). 

The pair functions on the Fermi surfaces ai^2 are 
isotropic and those on /3i^2 are anisotropic. The isotropic 
gap amplitude on ai_2 is Ao|$a|- The anisotropic gaps on 
/3i,2 have the maximum (minimum) value Aq (Ao$;3min)- 
From the RPA results presented in Ref. H^, it seems that 
^ 0.2 and ^fSmin ^ 0.2. With adjusting «'/3min, 
and Aq/Tc as parameters, we will calculate the nuclear 
magnetic relaxation rate 1/Ti and the superfluid den- 
sity pxx- We consider the following pair functions: (i) 
isotropic s-wave ($a < 0, <i>/3min = 1), (ii) anisotropic 
s-wave ($a < 0, ^pmin 7^ 1), (iii) isotropic ±s-wave 
(<i>a > 0, <i>/3min = 1), and (iv) anisotropic its-wave 
(<i>a > 0, <i>;3min 7^ 1)- Hcrc, wc cxcludc spiu-triplct pair- 
ings because Knight-shift measurements suggest a spin- 
singlet pairing42i^ 



III. NUCLEAR SPIN-LATTICE RELAXATION 
RATE 



The nuclear spin-lattice relaxation rate 1/TiT is given 



a* 



75.76.77.78 



(see Appendix) 



Ti{T)T AT 



duj 



cosh^(cj/2r) 



W{u;), (3.1) 



with 



Wiu;) = (aJ^(^)),3(a}i(-^)),,-(aJ|(^))^3(a|^(-..)),3. 

(3.2) 

= WGG{i^) + WFFiL0). (3.3) 
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Here, 



anisotropic s-wave 



ajj {kp, Uj) 



a]l{kp,uj) 



a'f\-{kF,u!) 



1 r 
2 



(fcp, ILUn 

g^l{kF,iuJn - 

-/n(feF,«w„ - 



f iT]) 
irj) 

irj) 
f «;) 
iri) 

Z7y) 



and 

5ll(feF,iu;„) 



.9ii(fcF,iw„) 
/l^(A;F,iw„) 
/xi(fcF,zw„) 



Vc.2 + |A(fcF)| = 

A(fcF) 



2 + |A(fcF)P 
A*(fcF) 
Vc.2 + |A(fcp)|2 



(3.4) 
(3.5) 
(3.6) 
(3.7) 

, (3.8) 
, (3.9) 
.(3.10) 



The brackets (• • • )fs rnean the Fermi-surface average, 

dSp i 



(...) 



E 

i=ai ,Q2,/3i./32 



FS 



E 

i=cti,a2,/3i./32 



d5] 



(3.11) 



k'F(feF) 



where dSp^i is the area elements on each Fermi surface. 
ujn = TTT{2n + 1) is the Matsubara frequency. We use 
units in which h — = 1- We assume the smearing 
factor T] = O.lTc.™ We set 2Ao/Tc = 4, which is a repre- 
sentative value near the BCS value 3.53. The coherence 
factor is represented as 1 + Wff/Wgg- The contribution 
of Wff is related to the coherence effect, which becomes 
zero in the case of unconventional pair functions such as 
d-wave one. 

First, we consider conventional s-wave pair functions 
(<I>Q < 0). In Fig. [21 we show the results for the isotropic 
s-wave ($0 — —1, *i'/3min = 1) and the anisotropic s- 
wave ($a = —1, *i'/3min = 0.2). The coherence peaks ap- 
pear below Tc for both pair functions because of non-zero 
Wff, meaning that these pair functions cannot explain 
the experiments. 

Second, we consider is-wave pair functions (<i>a > 0). 
Figure (Ha) shows the result in the case of the isotropic 
lbs-wave pair function (<i>Q — 1, $/3min = 1) and Fig.lD^b) 
shows the result in the case of the anisotropic its-wave 
function ($a = 0.2, $/3min = 0.2) whose fc-dependence 
is similar to the result of the RPA calculation by Arita 
et al^ In both cases, the coherence peak below Tc is 




FIG. 3: (Color online) Temperature dependence of the nu- 
clear magnetic relaxation rate 1/TiT (red circles) with the five 
band model in the case of (a) the isotropic s-wave ($„ = — 1, 
'I'/Jmin = 1) and (b) the anisotropic s-wave ($a = —1, 
^/Jmin = 0.2). 2Ao/rc — 4 and smearing factor 77 — O.lTc. 
The green squares denote the contribution from Wgg related 
to the density of the states and the blue triangles denote Wff 
related to the coherence effect. 



anisotropic ±s-wave 




FIG. 4: (Color online) Temperature dependence of the nu- 
clear magnetic relaxation rate 1 /T\T (red circles) with the five 
band model in the case of (a) the isotropic ±s-wave ($a = 1, 
'l'/3min = 1) and (b) the anisotropic ±s-wave ($a = 0.2, 
"l'/3min = 0.2). 2Ao/rc — 4 and smearing factor 77 — O.lTc. 
The green squares denote the contribution from Wgg related 
to the density of the states and the blue triangles denote Wff 
related to the coherence effect. 



suppressed, since Wff is almost zero. In the five band 
model, the difference of the density of states between the 
Fermi surfaces ai.2 and [3i^2 is small. Therefore, the can- 
cellation of the lbs-wave pair functions between a and (3 is 
almost perfect, resulting in Wff ~ 0. On the other hand, 
the temperature dependence at low T is inconsistent with 
the experiments in both cases. The exponential behavior 
appears in the isotropic ±s-wave case^i The tempera- 
ture dependence is concave down in the anisotropic ±s- 
wave case with $a = 0.2 and $/3inin = 0.2 as seen in Fig. 
UJb). We next consider another parameter set for the 
anisotropic is-wave pair function below. 

We search for the most suitable pair function with 
and $/3min- We check the two points as follows: (i) the 
lack of the coherence peak below and (ii) the low tem- 
perature power-law behavior 1/TiT oc T^. We show the 
temperature dependence of l/TiT in the cases of the var- 
ious pair functions in Fig. [S] First, we fix the /? gap 
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(a) anisotropic ±s-wave (b) anisotropic ±s-wave 




FIG. 5: (Color online) Temperature dependence of the nuclear 
magnetic relaxation rate 1/TiT with the five band model, (a) 
^0min = 0.2 and $a = 0.2 (1), 0.5 (2), 0.75 (3), and 1 (4). (b) 
= 1 and $^,nin =0.25 (5), 0.3 (6), 0.5 (7), 0.75 (8), and f 
(9). 2Ao/rc = 4 and smearing factor rj = O.lTc. The dashed 
line is a plot of T^. 



anisotropy <&^min — 0.2 and examine the (the a gap 
amplitude) dependence as shown in Fig. [5{a). With in- 
creasing $a from the value — 0.2, the exponent (i.e., 
the slope in Fig. [5]) approaches to the experimental result 
^ 2. The best coincidence is attained at $a = 1- Sec- 
ond, we fix <i>a = 1 and examine the ^/smin dependence 
as shown in Fig. [D^b) . With decreasing the anisotropy 
(i.e., increasing <I>/3min), the deviation becomes larger for 
*J?/3min > 0.25. Hcncc, the experimental results are best 
reproduced when = 1 and <i>^min = 0.25. That is, 
the maximum pair amplitudes on the Fermi surfaces ai^2 
and Pi^2 are of the same order ($a = 1), and the ratio 
of the minimum to the maximum of the pair amplitude 
on /3i_2 is 0.25 {^pmin — 0.25). We show the comparison 
of our calculation with the experimental result of ''''^As- 
NQR for LaFeAsOo.6 (Ref. M) in Fig. H Indeed, this 
anisotropic its-wave pair function explains the observed 
low-temperature power-law behavior l/Ti oc T^. 

IV. SUPERFLUID DENSITY 

Let us confirm whether the above anisotropic ±s-wave 
pair function can also explain the observed temperature 
dependence of the superfluid density. The superfluid den- 
sity pxx is given by^i^^ 

Pxx ^ 27rT / {t;Fx(fcF)}'|A(fcF)|" \ 

({v^k^)}') J^o \ K + |A(feF)r)'^' /fs' 

(4.1) 

Here, po denotes the superfluid density at the zero tem- 
perature and v-px is the Fermi velocity component in the 
(tt, 0) direction. 

As shown in Fig. [71 the superfluid density Pxx{T) 
for the anisotropic ±s-wave pair function (<i>£i 1, 
^/3min — 0.25) docs uot depend on the temperature in 
the low temperature region. When we increase 2Ao/Tc, 
the result approaches to that of the isotropic s-wave 
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FIG. 6: (Color online) Temperature dependence of the nu- 
clear magnetic relaxation rate l/Ti on a double- logarithmic 
scale. The red circles denote the result of the anisotropic its- 
wave pair function ($„ ~ 1, 'I'/gmin = 0.25, 2Ao/Tc — 4, and 
smearing factor r; — O.lTc). The green squares represent the 
experimental result of As-NQR for LaFeAsOo.e by Mukuda 
et al. (Ref. Hi). The dashed line is a plot of T^. Inset: Plots 
of the same data for 1/TiT on a non-logarithmic scale. 
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FIG. 7: (Color online) Temperature dependence of the su- 
perfluid density pxx for the anisotropic ±s-wave pair function 
(<E>a = 1, "J/jmin = 0.25). 2Ao/Tc = 4 (red circles), 5 (green 
squares). The dashed line represents pxx{T) for the isotropic 
s-wave gap with 2Ao/Tc = 4. po is the superfluid density at 
T = 0. 

case. Indeed, the anisotropic ±s-wave pair function 
can explain the fully-gapped behavior observed in the 
experimentsj ^^'^'^1^'^'^^1^^'^'' In contrast to it, pair func- 
tions with line nodes such as d-wave one lead to a strong 
temperature dependence near the zero temperature in 
general. 



V. DENSITY OF STATES 

Finally, we show the density of states N^{E) for the 
anisotropic is-wave pair function (<I>a = 1, $/3min = 
0.25) with 2Ao/rc = 4 in Fig. [8] It is calculated by 
N%E) = 7V"Re(5n(iw„ -> E where N"^ is the 
normal-state density of states at the Fermi level and g^^ 
is defined in Eq. p.Sp . The density of states is gapped in 
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FIG. 8: (Color online) Energy dependence of the density of 
states at T = for the anisotropic is-wave pair function 
(<E>a = 1, $^min = 0.25, and 2Ao/Tc = 4). A^" is the normal- 
state density of states at the Fermi level. Plots represent the 
density of states with the smearing factor rj — O.lTc (green 
dashed line) and O.OlTc (red sohd line).^^ 

the region \E\ < $/3minAo = 2<^>p„,inTc = 0.5Tc ($^„iinAo 
is the minimum gap on the Fermi surfaces /3i,2)- This is 
the reason why the superfluid density does not depend 
on the temperature in the low temperature region. In 
the region <^p^^^Ao{= 0.5T,) < \E\ < Ao(= 2T,), the 
density of states has a linear energy dependence. There- 
fore, the nuclear magnetic relaxation rate exhibits the 
line-nodes-like power-law behavior. The density of states 
also has the single peak structure near the gap edge at 
i?| = 2Tc = Aq, since the gap maxima on the Fermi sur- 
faces a 1.2 and /3i^2 now coincide with each other owing 
to $a = 1- Note here that the maximum gap amplitudes 
on ai,2 and /3i^2 are Ao|$a| and Aq, respectively. 

In addition, the density of states for the anisotropic 
lbs-wave pair function is a monotonically-increasing func- 
tion of the energy ([£'1 < Ag) as seen in Fig. [H while 
the unitary-scattering-induced density of states and the 
multi-gapped density of states are nonmonotonic in some 
cases. ^^'^^ This difference would be observed by spec- 
troscopy experiments. 



VI. CONCLUSION 

With the use of the five band model, we calculated the 
nuclear magnetic relaxation rate 1 /Ti and the superfluid 
density p^x and showed that the anisotropic is-wave pair 
function can explain the seemingly contradictory exper- 
imental results on Fe-pnictide superconductors. That is, 
the anisotropic is-wave pair function reproduces consis- 
tently 1/Ti ^ and the T- independence of p^x at low 
T. 

Our scenario is similar to the theories by Parker et 
al.^ Chubukov et al.,"^^ and Bang and Cho:'^^'""' in the 
sense that its-wave pair functions are considered in all 
theories. However, impurity effects are essential for those 
previous theories,^2i^i29 The impurity scattering rate is 
relatively large in Refs. W? and Issl . A unitary-limit im- 
purity scattering or an impurity scattering intermediate 



between Born and unitary limits^ is essential in Refs.lSTi 
and [3^ In contrast, we have assumed a rather clean sys- 
tem and not considered a unitary-limit or an intermediate 
phase-shift scattering. On the other hand, it was pointed 
out that a fitting resulted in quite big value 2Ao/Tc ~ 7.5 
within a model in Ref. [iol . In our model, rather strong 
gap anisotropy on the /? Fermi surfaces^ has been in- 
troduced, which enables us to explain 1/Ti even 
in a clean system and with relatively reasonable value 
2Ao/Tc ~ 4. This is a distinguished feature of our sce- 
nario. 

It should be noted that while some of experimen- 
tal groups have reported the fully-gapped behavior of 
the superfluid density, part of measurements showed 
somewhat strong temperature dependence indicating gap 
nodes^'^'2i!^i2^iS2.28,29,30 Those results seem to depend 

on kinds of materials and doping level, but it is still un- 
clear what is the essential origin of such scattered obser- 
vations between materials. The difference might mean 
that the pairing symmetry changes between materials or 
that the degree of gap anisotropy on the f3 Fermi surfaces 
changes, albeit there are no microscopic theories suggest- 
ing them at present. In any case, it is an interesting issue 
left for feature studies. 
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APPENDIX 

In this Appendix, we describe the procedure for de- 
riving the nuclear spin-lattice relaxation rate T^^{r,T) 
on the basis of the quasiclassical Green function 
theory^i^s^^^^^^ The derived formula has been 
utiHzed in Sec. Ill and in Refs. [zilTilT^lTi. 

Quasiclassical theory — We start with the Green func- 
tions defined as^ 

Gs^s'ir,r';T) = -(t, [V',(r, r)^^^ (r', 0)] ), (A.la) 

Fs,s' (r, r'; r) = -(t^ [^,(r, t)^;,, (r', 0)] ), (A.lb) 

F,,,,{r,r'-r) = -(t^VK^, t)V'1, (r', 0)] ), (A.lc) 

Gs,s'{r,r';T) = -(t^^J (^^ ^)^^, (^'^ q)] ). (A.ld) 
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Here, the brackets (■ • • ) denote the thermal average. We 
use units in which h = ~ 1. We write 




(A.2) 



Throughout this Appendix, "hat" {A) denotes the 2x2 
matrix in the spin space, and "check" (A) denotes the 
4x4 matrix composed of the 2x2 particle-hole space 
and the 2x2 spin one. 

The quasiclassical Green function g is defined as 



g = T3 



d^kG = T3 { Z21 ^22 



(A.3) 



where the integration is performed with respect to the 
energy variable in the k space, 



^A; = e{k) - fj.. 



(A.4) 



Here, s{k) is the quasiparticle dispersion relation and /i 
is the chemical potential. We have defined 



T3 



ao 
-<To 



with (To 



1 
1 



(A.5) 



According to a conventional procedure, the fc-space inte- 
gration is approximated as 

d^k ,r f f ,^ 
j^^N.J^Jd^k. (A.6) 

Here, an isotropic spherical Fermi surface is assumed 
for clarity. The extension to general cases can be done 
straightforward by replacing the solid-angle integration 
/ dil/ATT with the Fermi surface average (• ■ • )fs- is 
the total density of states at the Fermi level. 

The quasiclassical Green function follows the Eilen- 
berger equation, which is given as^^iS^iSiiS^iS^iSliS^iSa 

ivp ■ Vg + [iojnh - A, 5] = 0, (A. 7) 

where A is the superconducting order parameter, 

A = ( . " I . (A.8) 



A 

-At 



This equation is supplemented by the normalization 
conditional^ g'^ = -tt^I. 

We define, in the particle-hole space, the matrix ele- 
ments of the quasiclassical Green function g as^^ 



5„ 



(A.9) 



Comparing Eqs. (|A.3[) and (jA.9[) we have the following 
relation, which we will use later. 



11 



9 

9'' = 



9 



9'' 



nrg, 
ing, 
12 



(A.lOa) 
(A.lOb) 
(A.lOc) 
(A.lOd) 



In the case of spin-singlet superconductivity, the Eilen- 
berger equation is solved in a spatially uniform system 
and the solution for the quasiclassical Green function is2£ 



9 
f 



V^4 + W 

Ai&y 



V^4 + W 

V^^ + |AP' 



(A.ll) 



Here, the Pauli matrices are & = {ax, (jy^az) in the spin 
space. 

Relaxation Rate — The nuclear spin-lattice relaxation 
rate T^^(r,T) is obtained from the spin-spin correla- 
tion function x \-{x,x')fi^ We define x = {r,T), and set 

r' = 0. We apply a static external magnetic field along a 
certain axis and take the spin quantization axis parallel 
to this. X \-{x,x') is given as 

X-+{x,x') = (Tr[S-{x)S+{x')]) (A.12) 

= (Tr [^Pl (^)^T (^')^t i^')] j>A. 13) 
= Gii{x,x')G^^{x,x') 

~Fit{x,x')Fu{x,x'). (A.14) 

Let us consider a Fourier transformation with respect 
to T. In what follows, A and B stand for the Green 
functions. The Fermi- and Bose-Matsubara frequencies 
are ujn = ■nT{2n + 1) and f2„ — nT{2n), respectively. 
The Fourier transformation is 

A(r, r'; r) = ^ ^ M(r, r'; ic.„). (A.15) 



Note that A{t) and A{t)B{t) are periodic functions of 
r with the periods 2(3 and (3, respectively. Using Eq. 
()A.15p . we have the relation 

/ dTe'''^^A{T)B{T) = - VA(^c^„)B(^a„-^c^n). 

(A.16) 

From Eq. (|A.14p . the spin-spin correlation function is 

X-+(r,r';zr!„) = / dTe'f^'"-x-+(r, r'; r) (A.17) 
Jo 

1 

Gii{r,r';T)G^^{r,r';T) 

(A.18) 
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Using Eq. (fOBjl . we obtain 
X-+(r,r';ia^) = 4X1 

X ^Gii{r,r';iujn)G^^{r,r';inrn - iw„) 

(A.19) 

Now, we define f = r r' . 

*'=M(r,fe). (A.20) 



(A.21) 
(A.22) 



Setting iUm -^n + iS (5 ^0+), 



/OG nOO 
duj / dw'a^(w)a^(-cj') 
'OC J — OO 



. f{L0')-.f{L0) 



uj — uj' — — i5 
(A.28) 



duj 



aw a (tjja (— w ) — - 

UJ — uj' — It 



Air,r')^A{r,r)^ J ^^^^^ 
Setting r' ^ r (i.e., f = 0), we have 



Here, we have referred to Eq. (jA.6P (quasiclassical ap- 
proximation), k denotes the position of k on the Fermi 
surface. Then, from Eqs. (|A.19p and (|A.22|) . X-+(»': = 
r; iri,„) is 

X_+(r,r;ifJ„) 



(A.23) 



Here, we have referred to Eq. (jA.3p and have replaced 
/ d^/in with (• ■ • )fs- 

Next, let us consider the spectral representation of the 
quasiclassical Green functions: 



A{iuJn) 



- djio) 

OLU 



(A.24) 



Utilizing the formula (/(w) is the Fermi distribution hmc- 
tion) 



1 



we calculate 



(A.25) 



(A.26) 



duj 



doj'a'^{uj)a"{-Lo') 



(A.27) 



dio I duj' 



(A.29) 



1 



where we have used 
I 

Thus, 

Tmffl = . 

4 



(A.30) 



hm Imffl = ^ 

f2-»0+ VL 



duja^{Lj)a^{~uj) 



cosh^(/3t^/2) 



(A.31) 

It is known that Tf ^(r, T) is calculated by^^ {5 0+) 



TrVr,T) = T lim Ini 

0^0+ 



n 



(A.32) 

Referring to Eqs. C03)) . C06)) . COT]) , and (|X32)) . we 

obtain 

doj T 

, cosh^(w/2T) 



^ J— ( 



(«u(^))fs("Tt(-^))fs 
(«uH)fs(«u(-^)) 



(A.33) 

In the normal state, the spectral function of the quasi- 
classical Green function is a^^ — aP' = 1 for diagonal 
components in the particle- hole space (i.e., the density 
of states is unity in units of iVp) and is a^^ = a^^ = for 
off-diagonal components (because the order parameter is 
zero). We then obtain at T = Tc, 



duj- 



1 



cosh^(w/2rc) 



= TTTrN, 



(A.34) 
(A.35) 



Hence, the relaxation rate presented in Sec. HI is ob- 
tained: 



Ti{r,T)T AT 



duj- 



1 



cosh2(u;/2T) 
(«u(-))fs(«Tt(--))fs 



OC 

22/ 



(«uH)Fs(«n(— )) 



FS 



(A.36) 
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Spectral functions — In the spectral representation, the 
quasiclassical Green functions are 



duj 



IUJ„ — LO 



(A.37) 



where i, j = {1, 2}. Letting iw„ —>■ E ±iri (rj > 0), 



g^^ {iujn E ±iri) 



duj- 



E — u ±ir] 



(A.38) 



cLu 



a'J (lu) 
' E-uj 



TiTTa'^E), 



(A.39) 



where Eq. (jA.30|) has been used. From this, we have the 
relation 

Ztt l 

- g'^ {r,k,iu;n^ E-iT]) . (A.40) 
Referring to Eq. (jA.lOp . we have 



a^\r,k,E) = - 



a^^r,k,E) 



a^\r,k,E) = 



d^\r,k,E) 



g{r,k,iujn E + irj) 

g{r,k,iujn ^ E - IT]) , (A.41) 
g{r,k,iujn E + irj) 

g{r,k,iuOn E - irf) , (A.42) 
f{r,k,iuj„ E + iri) 

f{r,k,iun^ E -iri) , (A.43) 
/(r, k, icun E + irf) 

f{r,k,iun^ E -iri)]. (A.44) 



To calculate the relaxation rate in Eq. (|A.36[) . we need 
to consider the spin-space matrix elements presented in 
Sec. III. 
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